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Abstract. 

Novel  trapped  particle  asymmetry  modes  propagate  on  cylindrical  electron  columns  when  axial 
variations  in  the  wall  voltage  cause  particle  trapping.  These  modes  consist  of  E x B drifts  of  edge- 
trapped  particles,  partially  shielded  by  axial  flows  of  interior  untrapped  particles.  A simple  theory 
model  agrees  well  with  the  observed  frequencies  and  eigenfunctions;  but  the  strong  mode  damping 
is  as  yet  unexplained.  These  modes  may  be  important  in  coupling  trap  asymmetries  to  particle 
motions  and  low  frequency  E x B drift  modes. 


INTRODUCTION 

Single  species  plasmas  confined  in  cylindrical  traps  with  static  electric  and  magnetic 
fields  are  utilized  in  research  ranging  from  basic  plasma  and  fluid  dynamics  to  spec- 
troscopic frequency  standards  to  anti-matter  containment  [1].  The  plasma  confinement 
times  can  be  hours  or  days,  generally  limited  by  azimuthal  asymmetries  in  the  trap- 
ping fields  which  couple  angular  momentum  into  the  rotating  plasma  [2],  This  coupling 
can  be  enhanced  by  mode  resonances  [3],  and  plasma  modes  occur  on  frequency  scales 
ranging  from  the  (high)  cyclotron  frequency  £2C,  to  the  plasma  frequency  cop,  to  the  (low) 
E x B drift  frequency  (Of.  However,  the  well-studied  “diocotron”  drift  modes  have  axial 
wavenumber  kz  ~ 0,  and  so  do  not  readily  couple  the  plasma  electrostatic  energy  into 
particle  kinetic  energy. 

Here,  we  describe  novel  “trapped  particle  asymmetry”  modes  which  propagate  az- 
imuthally  at  low  E x B drift  frequencies,  but  which  have  kz  0 and  thus  couple  to 
particle  kinetics.  These  modes  exist  when  axial  variations  in  the  wall  potential  cause  the 
equilibrium  plasma  to  have  axially  trapped  particles.  Experimentally,  we  generate  the 
trapped  particles  by  applying  a “squeeze”  voltage  to  a central  cylindrical  electrode;  but 
smaller  trapping  potentials  are  caused  by  wall  irregularities,  and  are  probably  endemic 
to  these  traps. 

These  asymmetry  modes  can  have  azimuthal  mode  number  m — 1,2,...,  but  we 
focus  on  m = 1 here.  Dynamically,  the  modes  are  an  E x B drift  motion  of  the  edge- 
trapped  particles,  with  axial  flows  of  interior  untrapped  particles  giving  partial  Debye 
shielding.  A simple  theory  model  with  these  characteristics  shows  close  correspondence 
with  the  measured  frequencies  and  radial  eigenfunctions  of  the  modes.  The  modes  are 
observed  to  be  strongly  damped,  but  this  damping  is  not  yet  understood.  These  modes  are 
particularly  simple  cylindrical  analogues  to  the  trapped  drift  modes  which  can  contribute 
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FIGURE  1.  Schematic  of  the  cylindrical  trap  with  "squeezed”  electron  column  and  applied  sector 
voltages. 


to  anomalous  plasma  transport  in  tokamak-like  and  toroidal  multipole  configurations  of 
neutral  plasmas  [4,  5], 


EXPERIMENTAL  SETUP 

The  cylindrical  confinement  geometry  and  a schematic  of  the  mode  are  shown  in  Fig.  1. 
A nominally  cylindrical  column  of  electrons  from  a hot  tungsten  source  is  confined 
radially  by  a uniform  magnetic  field  B = 4 kG  and  confined  axially  by  negative  voltages 
-Vr  = -100  V on  end  cylinders.  Typical  electron  columns  have  a central  density  no  « 
1.5  x 107  cm"3  over  a length  Lp  « 50  cm,  with  a column  radius  of  Rp  w 1 .2  cm  inside  a 
wall  radius  Rw  = 3.5  cm.  The  z-averaged  electron  density  «(r,0,r)  can  be  (destructively) 
measured  at  any  time  by  dumping  the  electrons  axially  onto  a phosphor  screen  imaged 
by  a CCD  camera  [6], 

The  individual  electrons  have  a thermal  energy  T ss  1 eV,  giving  an  axial  bounce 
frequency  = v/2Lp  « 0.5  MHz  and  a Debye  shielding  length  Xo  = ( T /4ne2n)^2  « 
0.2  cm.  The  unneutralized  electron  charge  results  in  a central  potential  <J>o  « -30  V, 
and  the  radial  electric  field  causes  the  column  to  E x B drift  rotate  at  a rate  /£■(/■)  = 
cE(r)/2nrB  < 50  kHz.  Figure  2 shows  typical  radial  profiles  of  density  and  rotation  for 
a 0-symmetric  column. 

When  a negative  “squeeze”  voltage  -Kq  is  applied  to  a central  cylinder,  electrons  are 
excluded  from  the  column  periphery  under  the  squeeze  ring,  and  those  electrons  located 
at  radii  at  r > rs  are  trapped  axially  in  one  end  or  the  other.  For  small  VSC|,  a fraction 
Etr  ~ Esq/tyo  l*ie  electrons  are  trapped  axially. 

Linear  E x B drift  modes  with  azimuthal  mode  numbers  m—  1,2, ...  propagate  on  this 
trapped  particle  equilibrium.  The  ubiquitous  m = 1 center-of-mass  “diocotron”  mode 
is  essentially  uniform  in  z (i.e.  kz  & 0),  and  is  nominally  unaffected  by  small  Fsq.  In 
contrast,  the  m = 1 “trapped  particle  asymmetry”  mode  described  here  has  odd  parity 
in  z:  the  perturbations  are  essentially  uniform  within  each  end,  but  of  opposite  sign  on 
either  side  of  the  squeeze,  as  shown  schematically  in  Fig.  1 . Moreover,  the  perturbations 
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FIGURE  2.  Radial  profiles  of  density  n(r),  and  E x B rotation  frequency  /Hr). 


in  the  trapped  particles  at  r > rs  are  partially  “shielded”  by  perturbations  of  opposite 
sign  in  the  untrapped  particles  at  r < rs\  and  these  untrapped  particles  slosh  from  end  to 
end  in  response  to  the  E x B drift  evolution  of  the  trapped  particles. 

Experimentally,  these  odd-z-parity  m = 1 asymmetry  modes  can  be  excited  with  + 
and  - voltages  applied  to  opposing  wall  sectors,  with  the  voltages  oppositely  phased  at 
the  two  ends.  Any  single  wall  sector  can  be  used  as  a receiver,  or  sectors  can  be  used  in 
combination  to  verify  or  utilize  the  0-  and  z-symmetries  of  the  modes. 


EXPERIMENTAL  RESULTS 

Figure  3 shows  the  measured  frequencies  /«  and  fj  as  the  applied  squeeze  voltage  Esq 
is  varied.  The  asymmetry  mode  frequency  fa  is  at  or  near  the  edge  rotation  frequency 
/e(^p)  for  ESq  -C  <t>o-  As  Esq  is  increased,  the  trapping  separatrix  rs  moves  inward,  and 
fa  decreases.  For  Esq  > <j>o,  the  column  is  cut  in  half,  and  the  asymmetry  mode  becomes 
degenerate  with  the  diocotron  mode.  With  all  particles  trapped,  the  asymmetry  mode  is 
equivalent  to  two  separate  diocotron  modes,  180°  out  of  phase.  The  diocotron  frequency 
fd  increases  slightly  with  VSq,  because  the  effective  line  density  of  the  column  increases 
as  particles  are  excluded  from  the  squeeze  ring.  Both  fa  and  fa  scale  as  /°=  B~l,  as 
expected  for  E x B drift  modes. 

Figure  4 shows  that  the  asymmetry  mode  damping  rate  ya  also  depends  on  Fsq, 
decreasing  roughly  as  ya  (1  — Esq/<(>o)  in  the  regime  Vsq  > <j)o/10  where  ya  can  be 
readily  measured.  Comparing  Figs.  3 and  4,  we  can  see  that  y a/ fa  1 /20,  so  the  mode 

is  strongly  damped  except  at  Esq  > <|)o.  However,  we  observe  essentially  no  magnetic 
field  dependence  to  the  damping  [at  fixed  Esq/<j>o]  in  the  range  of  0.39  kG  < B < 10  kG, 
i.e.  ya  B°,  (A  weak  dependence  on  B arises  experimentally  because  the  equilibrium 
plasma  potential  <t>o  depends  weakly  on  B.)  Thus,  y a/fa  x B]  and  the  mode  is  relatively 
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FIGURE  3.  Measured  frequencies  /„  and  /,/  vs  squeeze  voltage  VV  fa  calculated  from  the  kinetic 
model  (dashed). 


Ya  [sec1] 


FIGURE  4.  Measured  decay  rate  y,i  of  the  asymmetry  mode  vs  squeeze  voltage  FSq  at  two  extreme 
values  of  the  magnetic  field.  Plasma  potential  <{>o  is  less  at  low  magnetic  field  due  to  enhanced  radial 
transport  of  the  particles. 
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FIGURE  5.  Measured  (circles)  and  calculated  from  the  kinetic  model  (solid)  asymmetry  mode  eigen- 
functions 8 na(r)  for  VSq  = — 5V  (left)  and  Vsq  = — 36V  (right). 


less  strongly  damped  at  lower  magnetic  fields. 

As  yet,  the  cause  of  this  damping  is  not  understood,  but  experimental  signatures  (like 
ya  rjfi0)  suggest  diffusive  mixing  between  the  trapped  and  untrapped  populations. 
Also,  for  Vsq  < <J)o/10,  there  are  suggestions  that  spatial  Landau  damping  (scaling  as 
fi_1)  may  occur,  but  this  damping  depends  critically  on  the  density  profile  near  Rp. 

The  density  eigenfunctions  8 na(r)  of  the  asymmetry  mode  can  be  obtained  from  a 
time  sequence  of  measurements  of  the  density  n/,(r,6,t)  z- averaged  over  only  half  the 
plasma.  To  dump  only  half  the  column,  Vsq  is  increased  to  -100  V immediately  before 
lowering  Vc  to  dump  the  plasma  onto  the  phosphor  screen.  The  m = 1 perturbations  in 
«/,(/-,  0,/)  are  then  fit  to  a sum  of  2 modes  (asymmetry  and  diocotron),  as 

[ dd  nh{r,  0,  t)  e~imQ  = £ &i;(r)  e'2n//  ~ Y j* . 

j—a.d 

The  residual  to  this  fit  is  small,  and  the  eigenfunctions  are  obtained  with  little  ambiguity. 
The  real  parts  of  the  asymmetry  mode  eigenfunction  ?ma  (scaled  arbitrarily)  obtained 
for  — — 5V  and  Vsq  = —36V  are  shown  in  Fig.  5;  with  proper  choice  of  0-origin,  the 
imaginary  part  of  8 na  is  essentially  zero. 

The  mode  shows  a negative  perturbation  at  r > rs  and  a positive  perturbation  for  r <rs\ 
and  the  signs  of  these  perturbations  reverse  for  z — > — z and  for  0 — > 0 + n.  The  observed 
position  of  the  zero-crossing  at  rs  varies  with  Vsq,  consistent  with  the  “waistline”  radius 
obtained  from  (r,z)  Poisson  solutions  for  the  equilibrium  with  the  measured  charge  and 
applied  wall  potentials. 
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THEORETICAL  MODELS 


Fluid  theory  and  kinetic  theory  models  further  support  the  inteipretation  of  the  asymme- 
try mode  as  trapped  E x B drifting  edge  particles  which  are  Debye-shielded  by  sloshing 
interior  particles.  In  these  models,  the  applied  squeeze  voltage  is  presumed  to  create  a 
(zero-length)  barrier  for  particles  at  r > rs,  creating  two  outer  regions  with  axial  lengths 
L\  and  Z/>.  In  contrast,  interior  particles  with  r < rs  move  freely  over  length  L\  + Lt- 
Because  the  axial  bounce  frequency  is  large  compared  to  the  E x B rotation  frequency 
and  the  mode  frequency,  the  bounce-average  density  perturbation  in  any  one  of  the  three 
regions,  5/;,  is  related  to  the  bounce-average  potential  perturbation  in  that  region,  64), 
through 


Mr) 


c df/Q  »i8<t>(r) 
2v.Br  dr  mfE(r)-f  ' 


(1) 


Within  any  of  the  regions,  rapid  axial  streaming  yields  the  adiabatic  response 


8n{r,z)  = 8n(r)  + — [ 8<|>(r,z)  -8<|>(r)  ] • 


(2) 


This  type  of  response  gives  rise  to  Debye  shielding,  making  S<|>(r, z)  and  8n(r,z) 
nearly  z-independent  except  near  the  plasma  ends  and  the  squeeze  region.  Thus,  we 
approximate  8 n(r,z)  « 8 iij(r)  and  8<})(r,z)  8cJ)7 (r),  with  j — 1,2  representing  the  left 

and  right  sides  relative  to  the  barrier.  Poisson’s  equation  for  the  two  sides  is  then 


1 a 

r dr' 


2 ec  diiQ  m8§j 
Br  dr  mfE  - f 

2ec  3/io  /«8(J)/,  4ne2n  , 


,r>rs 
, r < r. 


(3) 


Here,  8(J)*  = (Ti8([>i  ^284)2 )/(Z-i  + Lj)  is  the  potential  averaged  over  both  sides,  and 

the  boundary  condition  is  84 >;(/?„  ) = 0. 

Fortunately,  the  even  and  odd  axial  symmetries  of  the  diocotron  and  asymmetry 
modes  decouple  the  equations.  For  84>i  {>')  = 54>2('  ) = 8<t>t(r),  both  equations  reduce  to 
the  usual  equation  for  a diocotron  mode.  For  Lj  84>i  (r)  = — IsSMr),  giving  8 <(> /,(/•)  = 0, 
we  obtain  an  eigenvalue  equation  for  the  asymmetry  mode  potentials. 

For  a uniform  density  profile  of  density  hq  and  radius  Rp,  the  eigenfrequency  / = f, 

fa  + 

/£(0)  m 

Here,  l,„±\  = lm±\{rs/Xo)  are  modified  Bessel  functions  of  the  first  kind.  In  the  limit 
of  Xp  ~ > 0,  the  interior  region  acts  like  a conductor,  and  the  exponentially  large  Im±\ 
may  be  replaced  by  1 in  Eq.  (4).  The  boundary  rs  is  determined  by  radial  integration 
of  Poisson’s  equation  with  density  n 0 out  to  radius  ry  giving  <| )(/?„,)  = — VS(|.  Figure  6 
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FIGURE  6.  Separatrix  radius  rs  vs  Vsq/<t>o  for  uniform  density  model,  giving  /„  from  Eq.  (4)  and 
constant  fa. 

shows  rs,  fa,  and  fa  versus  Vsq  for  this  uniform  density  model  with  Rp  — 0.5 RiV  and 
m = 1.  There  is  an  obvious  general  correspondence  between  this  coarse  estimate  and  the 
experiment  (Fig.  3). 

A kinetic  analysis  allowing  for  trapped  and  untrapped  particles  co-existing  at 
each  radius  gives  a more  realistic  approximation  to  the  experiments.  The  squeeze 
causes  a potential  barrier  of  strength  A<t>(r),  giving  a trapped  particle  density 
n,(r)  — n(r)erf[(<?A<t)/7')1/2]  and  an  untrapped  density  nu  = n — n,.  The  equilibrium 
Poisson  solutions  show  that  A<[>  is  essentially  zero  for  small  radii  and  that  it  rises  sharply 
near  rs  to  a value  much  larger  than  T.  Thus,  the  kinetic  treatment  essentially  smoothes 
the  discrete  transition  model  over  a radial  scale  of  a few  XD. 

The  frequency  fa  versus  Vsq  predicted  by  this  kinetic  theory  using  the  measured  n(r) 
and  the  calculated  A<j>(r)  is  given  by  the  dashed  curve  in  Fig.  3,  showing  close  agreement 
with  the  measured  frequencies.  The  10%  discrepancy  may  reflect  the  non-zero  length  of 
the  barrier  and  consequent  increase  in  trapped  particle  densities. 

Similarly,  the  eigenfunctions  5 n„(r)  obtained  from  the  kinetic  analysis  show  general 
correspondence  with  experiments  (Fig.  5).  The  discrepancy  for  r > 1.5  cm  may  be  an 
artifact  of  the  profile  truncation  applied  theoretically. 


MAGNETIC  TILT,  MODE  COUPLING  AND  TRANSPORT 

Although  the  asymmetry  modes  can  propagate  in  a trap  with  perfect  9-symmetry,  our 
larger  interest  is  in  traps  with  (inevitable)  0-asymmetries.  These  asymmetries  commonly 
arise  from  construction  imperfections  in  the  trap  walls  (especially  in  the  sectored  cylin- 
ders) or  in  the  magnet.  The  asymmetries  can  have  arbitrary  0-  and  z-dependence,  but  it 
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FIGURE  7.  Measured  frequency  /„  and  decay  rate  y,,  vs  magnetic  tilt  0g,  normalized  to  values  at  0«  = 0. 


is  common  to  refer  to  their  dominant  m and  k-  components.  One  of  the  more  prevalent 
asymmetries  is  a in  — 1 , kz~  n/Lp  “lilt”  of  the  magnetic  field  with  respect  to  the  trap 
axis,  characterized  by  an  angle  0#  = BA/B:. 

The  qualitative  behavior  of  the  asymmetry  mode  remains  unchanged  even  with  mod- 
erately large  magnetic  tilt.  Figure  7 shows  /„  decreasing  linearly  by  30%  and  y„  decreas- 
ing by  60%  at  fixed  Vsq  = -10  V as  the  magnetic  tilt  is  increased  to  0#  = ±3  x 10  3 
radian.  These  quantitative  changes  are  mainly  due  to  a decrease  in  the  plasma  potential 
<J>0  (caused  by  tilt-induced  particle  transport).  Consequently,  the  ratio  VSq/<},()  increases 
and  the  effective  trapping  barrier  is  enhanced.  This  will  produce  a decrease  in  the  mode 
frequency  and  decay  rate  according  to  the  trends  depicted  in  Figs.  3 and  4. 

Trapped  particles  may  be  an  inherent  part  of  Penning-Malmberg  traps  for  a variety  of 
reasons.  Some  intentional  manipulations,  such  as  “squeeze  damping”  of  the  diocotron 
mode  [7,  8]  obviously  create  equilibria  as  described  here.  The  wide  variety  of  “nested 
traps”  being  built  for  overlapping  confinement  of  positrons  and  anti-protons  [9,  10]  have 
populations  of  trapped  and  untrapped  particles.  Most  subtly,  z-variations  in  the  effective 
wall  voltage  can  arise  from  small  unintentional  construction  anomalies.  For  example,  a 
variation  in  wall  radius  among  cylinders  gives  an  effective  potential  v/sq  « <>„  A Rw/Rn.- 
this  gives  a trapped  fraction  e,r « A Rw/Rw,  so  8ir  > 10  3 is  probably  common  to  all 
devices. 

The  existence  of  trapped  particle  asymmetry  modes  substantially  alters  the  theory 
perspective  on  particle  transport  due  to  0-asymmetries  in  the  trap  construction.  Reso- 
nant particle  transport  theories  [3]  utilize  integration  along  unperturbed  orbits,  which 
is  invalid  with  the  bulk  trapped  particle  populations  described  here,  or  even  with  mi- 
croscopic trapping  [11],  Similarly,  2D  bounce-averaged  invariants  suggest  that  particles 
are  confined  to  equipotential  surfaces,  and  a kz  variation  is  required  for  radial  particle 
transport;  however,  all  previously-known  modes  with  kz  ^ 0 have  frequencies  near  co/; 
as  opposed  to  fj. 

Experimentally,  there  are  strong  indications  that  the  asymmetry  mode  contributes  to 
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FIGURE  8.  (a)  “Excitability”  of  the  asymmetry  (closed  circles)  and  diocotron  (open  circles)  modes  vs 

squeeze  voltage  Vsq.  (b)  Measured  ratio  fa/fd  (squares)  vs  squeeze  voltage  vv  The  horizontal  dashed 
lines  show  the  resonance  [Eq.  (6)]  of  sloshing  particles.  The  vertical  dashed  line  indicates  an  obvious 
correspondence  of  these  resonances  to  the  maxima  in  the  diocotron  mode  amplitude. 


asymmetry-induced  transport  by  coupling  between  the  static  m / 0 kz^0  asymmetry 
and  the  kz  0 diocotron  modes.  When  the  external  asymmetry  is  a magnetic  tilt  (0#  ^ 
0),  the  diocotron  mode  has  axially  sloshing  particles  [12]  at  low  frequency  fd  -C  fb, 
which  can  be  thought  of  as  a slow  change  in  the  equilibrium.  A simple  picture  of 
what  occurs  with  a misaligned  plasma  column  is  that  the  particles  tend  to  accumulate 
wherever  they  can  be  closer  along  the  magnetic  field  lines  to  their  image  charge  in 
the  wall,  since  that  is  a region  of  lower  energy.  Consequently,  a small  amount  of  the 
density  (8 n °c  0B)  sloshes  axially  from  one  end  to  the  other  as  the  diocotron  mode 
causes  the  plasma  column  to  orbit  around  the  magnetic  axis  of  the  trap.  This  axial 
motion  can  then  couple  to  the  axially  sloshing  interior  particles  in  the  asymmetry  mode. 
The  coupling  leads  to  the  diocotron  mode  damping  through  excitation  of  the  strongly 
damped  asymmetry  modes,  and  consequently  to  bulk  radial  transport  of  the  particles. 

The  characteristics  of  this  coupling  are  presently  being  investigated.  Figure  8(a)  shows 
the  “excitability”  of  the  asymmetry  and  diocotron  modes  as  functions  of  the  applied 
squeeze  voltage.  “Excitability”  is  simply  the  initial  amplitude  of  the  mode  excited  from 
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the  wall  sectors  by  a fixed  m = 1 kick  with  odd  parity  in  z (Fig.  1).  The  small  excitation 
of  the  diocotron  mode  at  Fsq«0  occurs  only  because  the  sectors  at  ±z  have  slightly 
different  coupling  to  the  plasma.  As  squeeze  voltage  is  increased,  the  asymmetry  mode 
“excitability”  increases  linearly,  in  accord  with  the  change  of  the  trapped  particle  fraction 
etr.  In  contrast,  the  excitation  of  the  diocotron  mode  shows  increasing  oscillations  about 
the  Vsq  w 0 level  [Fig.  6(a)], 

The  maxima  in  these  oscillations  may  be  due  to  resonance  in  the  axial  velocities  of 
sloshing  particles  shared  by  the  two  modes.  Thus,  we  could  expect  that  in  resonance 

/e  - fa  _ ft:  - tfd 

2nka  2nkd  ’ ' 

where  ka  and  kd  are  effecctive  axial  wave  vectors  ( ka  > kd)  for  the  sloshing  particles  in 
corresponding  modes,  and  (.  = 2,3...  is  the  order  of  the  resonance.  Figure  8(b)  shows  the 
measured  ratio  /„/ fd  as  a function  of  the  squeeze  voltage.  Comparing  Figs.  8(a,b),  we 
see  that  the  maxima  in  the  diocotron  mode  amplitude  corresponds  to  the  particular  ratios 
of  these  mode  frequencies,  obtained  directly  from  Eq.  (5)  using  a quite  reasonable  ratio 
of  the  wave  vectors  kajkd  « 1.45  as  the  single  adjustable  parameter 

k = & + &.(  1-M.  (6) 

Jd  «</  J,i  \ klt ) 

This  result  may  provide  a clue  as  to  why  the  diocotron  mode  has  shown  a squeeze 
dependent  damping  [7,  8],  Of  course,  the  underlying  assumptions  and  simplifications 
of  the  observed  coupling  remain  to  be  examined.  Nevertheless,  our  results  allow  us 
to  speculate  that  the  asymmetry  modes  are  important  in  coupling  trap  asymmetries  to 
particle  motions  and  low  frequency  E x B drift  modes. 
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